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Abstract 

We perform a calculation of the isospin and hypcrcharge vector current 
propagators {Ayi^^lq^} and Aygg(g^)) to two loops in chiral perturba- 
tion theory. The analysis is carried out with straightforward Feyn- 
man diagram methods by making appropriate use of external vector 
sources. Counterterms from the 0{q^) chiral lagrangian, required to 
absorb divergences and scale dependence encountered at the two-loop 
level, are constructed. Our final results are finite, covariant, and scale- 
independent. Several applications are described, including a compari- 
son of the two-loop isospin vector spectral function with data and the 
construction of new chiral sum rules. 
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1 Introduction 



A decade ago, Gasser and Leutwler performed a general one-loop analysis 
of chiral perturbation theory. ffl'B Among their results were representations 
of the isospin vector current and axialvector current propagators, valid up to 
0{q'^) corrections. In this paper, we shall describe a calculation to two-loop 
order in chiral perturbation theory^ (hereafter ChPT) of both the isospin 
and hypercharge vector current propagators, yielding representations valid 
up to 0{q^) corrections. 

In the years following the publication of the Gasser-Leutwyler papers 
cited above, there appeared a dramatic increase in the number of stud- 
ies utilizing ChPT to analyze hadrons at low energies. However, most of 
these involved either tree-level or one-loop calculations. In fact, there cur- 
rently exists a need for higher-order calculations which meaningfully probe 
the convergence properties of ChPT. This forms an important impetus for 
performing our calculation. 

There is a more phenomenologically-oriented way of stating our motiva- 
tion which actually bears on future applications of our work - the subject 
of chiral sum rules. It has long been known that among the predictions of 
chiral symmetry are the following set of sum rules, 

oo 

ds {py{s) - pa{s)) = , (2) 

oo 

ds s {py{s) - pAis)) = , (3) 
^ ds slni-^j {py{s) - pa{s)) = ^^{ml± -mlo) . (4) 

In the first of thesefl, the quantity L]|q(M^) is the renormalized coefficient 

of an 0{E'^) counterterm operator appearing in the SU{3) chiral lagrangian 
of QC-D.B We do not comment further on this quantity here insofar as 
much of Sect. 4 and Sect. 5 is devoted to a study of counterterms and their 
renormalization. Thenext two relations are respectively the first and second 
Weinberg sum rules. O The final sum rule is a formula for the vr^-vr'^ mass 

■^We stress for clarity's sake that our analysis throughout will be in the standard ChPT 
and not the so-called generalized ChPT of Ref. B . 
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splitting in the chiral limit. B Although apparently containing an arbitrary 
energy scale A, this sum rule is actually independent of A by virtue of the 
preceding relation. 

Although the above sum rules were derived almost 30 years ago, a full 
phenomenological test of their content was not carried out until recently^ . 
In principle, the underlying strategy employed in Ref. was to divide the 
energy region into three intervals as follows. In the large energy region, 
from s ^ oo down to energies where data first become available, operator 
product expansion methods were used. At intermediate energies, the vector 
and axialvector spectral functions py{s) and pa{s) were extracted from tau 
decay distributions and e~^e~ — > 27r, Att cross sections. Finally, the threshold 
behavior (s — > 4m^) of the spectral functions was constrained by means of 
chiral symmetry. In practice, however, this final step was not found to be 
crucial to the calculation and as a result was not implemented. The analysis 
of Ref. ultimately led to a successful fit to all four chiral sum rules, 
which is highly nontrivial in view of the strong distortions experienced by 
the spectral functions from the various moments. 

Let us study the threshold prediction of ChPT for the isospin spectral 
function. The dashed curve in Fig. 1 is the one-loop prediction of Gasser 
and Leutwyler and is to be compared with data taken from e^e~ tt^tt^ 
scattering. One readily sees that it does not take long for the one-loop 
prediction to disagree with experiment as the energy is increased above the 
two-pion threshold. This is one of several reasons for computing to the next 
order of ChPT.tEl 

We conclude this section with an outline of the rest of the paper. Section 
2 contains a discussion of our overall strategy for carrying out the calcula- 
tion. The lowest-order lagrangian C^'^^ is used to illustrate certain aspects 
of the discussion. Section 3 summarizes both the one-loop and two-loop 
chiral analysis and introduces the 0{q^) chiral lagrangian, to the cal- 
culation. Section 4 covers the counterterm analysis, reviewing the presence 
of the 0{q'^) counterterms and especially generating a complete list of 0{q^) 
counterterms relevant for our purposes. Section 5 contains a statement of 
our renormalization procedure. The use of the renormalization group (RG 
hereafter) to determine the scale dependence of various counterterms is dis- 
cussed and the RG equations are solved to the needed order. Section 6 
contains the explicit representations for the isospin and hypercharge vector 
current propagators which result from our calculation. Section 7 contains 
two applications of our results. The improved isospin vector spectral func- 



2 



tion, now containing next-to-leading order contributions, is again compared 
with data. In addition, the asymptotic behaviour of the vector-propagators 
is employed to derive a new set of chiral sum rules. We conclude the main 
body of the paper in Section 8 by summarizing our findings and proposing 
further research on the subject. There are also two Appendices, the first 
devoted to defining the d-dimensional integrals which occur in our calcula- 
tion and the second to giving an independent determination of the isospin 
spectral function py33{s) at next-to-leading order. 

2 Procedure 

Consider the SU{3) vector currents, @ 

Va^ = q^rq , (5) 
where a = 1, ... ,8 and q = [u d s). For the cases a = 6 = 3 and a = b 



[10 



the corresponding two-current time-ordered products can be expressed as 
{0\TiV,^{x)V,^iO))\0) = 

roc f ^4 p-ip-x 







where the quantities p\ah{s) are the spin-one vector spectral functions. The 
vector current propagators are Fourier transforms of the spatial two-point 
functions, 

KIM^) fd'x e'''-^{0\T{V,^{x)Vi'{Om . (7) 



In this paper, we shall use ChPT to calculate the low energy behavior 
of the isospin and hypercharge vector-current propagators Ay'^^{q'^) and 
Aygg(g^). Since the underlying current is conserved for each of these, we 
have the tensor decomposition 

Ara6(9') = ('?V-gV")na6(g') (a = 6 = 3, 8) . (8) 

We shall refer to this tensor structure as covariant. In the course of calculat- 
ing Ay33(g^) and Ay-gg(g^), one encounters a different class of contributions 
which are proportional to g'^'^ . We call these terms noncovariant. Of course, 
the set of noncovariant terms must ultimately cancel. 
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External Vector Sources 



In accordance with the plan to carry out a general analysis of ChPT 
to one-loop, Gasser and Leutwyler employed the so-called background-field 
method. yJ This had the advantage of generating all possible 0{q^) countert- 
erms. However, since our attention will be strictly limited to the conserved 
vector-current propagators U^s and IIss, we have adopted a rather different 
approach. We observe that it is sufficient to compute all Feynman diagrams 
of the type depicted in Fig. 2, i.e. those with a single incoming external 
vector source undergoing a transition to a single outgoing external vector 
source. It will become clear that this method is both accessible and straight- 
forward to implement. 

Let us illustrate the procedure in terms of the lowest-order chiral la- 
grangian 

£(2) = ^ Tr (d.UD'^U^) + ^ Tr [xU^ + Ux^) . (9) 

In the above, Fq is the pseudoscalar meson decay constant to lowest order 
and X = '^Bqui is proportional to the quark mass matrix with 



Bo = ^ = ^^= , . (10) 



Note that we work in the isospin symmetric limit of = = m. The 
field variable U is defined in terms of the pseudoscalar meson fields {(j)a}, 

U = exp(iAa • ^a/Fo) , (11) 

and we construct the covariant derivative D^U in terms of external vector 
sources u^, 

D^U = d^U + iv^U - iUv^ . (12) 
The vector source has a component v^^ for each of the SU (3) flavors, 

v,^Y'< ■ ^^^^ 

Suppose that we use the above form of C^'^^ to compute a Feynman 
diagram of the type in Fig. 2. For now, consider a term in C^"^^ which is 
linear in the vector source v^, 

Ant = -<(x)Knx) • (14) 
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The quantity V^{x) which is coupled to the external source v'^ must, by 
construction, be the corresponding vector current. To calculate the given 
Feynman diagram, we must work with the second-order operator ^2, 

S2-l = jJ<fx (fy T (£int(x)/:int(y)) . (15) 

Upon taking the matrix element for the diagram in Fig. 2 and employing 
translation invariance, we easily obtain the result, 

-l\i)= z(2^)45W(g' - q) el{q')MZeM , (16) 
where the invariant amplitude is found to be 

MZ = i J d'xe'^-^0\T{V,^{x)V,^iOm ^ K^abi^') ■ (17) 

That is, one just computes the invariant amplitudes for the class of diagrams 
in Fig. 2 and the result is the vector-current propagator of the relevant 
flavor. The only complication is that amplitudes corresponding to the linear 
interaction of Eq. (|T^ ) are generally noncovariant and one must thus include 
contact terms which are bilinear in the external vector sources. For example, 
in SU{3) notation the full interaction induced by the chiral lagrangian 
is 

Cint = -fabc(^bd^<l)cVl + -fadgfbcgi^cCpdV^v'^ , (18) 

where the two terms are to be calculated respectively in second and first 
order perturbation theory. In the next section, we shall use this lowest- 
order interaction to begin our treatment of the vector-current propagators 
to one-loop order. 

To recap briefly, our calculational program is to use isospin and hy- 
percharge vector sources and Vg as aids in computing the corresponding 
vector-current propagators. We shall require the chiral lagrangians 

£(2)^ £(4) 

and £^^\ These lagrangians will involve just the so-called normal vector cur- 
rents. In this paper, we will not consider any of the anomalous currents. In 
the course of computing the one-loop and two-loop Feynman diagrams, we 
shall encounter divergent integrals to which we shall assign precise meaning 
via dimensional regularization. Due to the many terms which occur, the 
calculation is both difficult and complicated. Fortunately, it will be possible 
to check the correctness of our results along the way in a variety of ways. 
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3 One-loop and Two-loop Analyses 



We shall divide discussion of our calculation into three separate parts — loop 
analysis, counterterm analysis, and renormalization procedure. This section 
begins the process with a description of one-loop and two-loop amplitudes 
which are obtained by using the vector-source approach of Section 2. The 
Feynman diagrams underlying our loop and counterterm analyses are the 
ones given in Fig. (3) and Fig. (4). Throughout, we shall denote all contribu- 
tions from the chiral by a small circle enclosing a cross and those from 
by a small square enclosing a cross. It turns out that all the Feynman 
integrals we shall encounter can ultimately be expressed in terms of just 
two quantities, which we call A and B21. They are defined in Appendix A, 
where we also discuss in detail all such integrals occurring in our analysis. In 
our dimensional regularization treatment of these quantities, we have been 
very careful to allow throughout for sufficiently high powers in the small- 
ness parameter d — 4. As it happens, the only terms in our expansions which 
actually contribute are the ones we display in App. A. 



The one-loop chiral analysis of the isospin vector-current propagator 
was originally carried out by Gasser and Leutwyler in an SU{2) basis of 



fields^. We shall briefly summarize the theoretical analysis with the SU{3) 
basis employed in this paper and shall also include the one-loop expression 
for the hypercharge propagator. Our discussion of the one- loop sector will 
serve to set the stage for the two-loop analysis. 

At one-loop level, the vector-current propagators are determined from 
the Feynman diagrams appearing in Fig. 3. There are two amplitudes, each 
arising from the chiral lagrangian C^'^\ the unitarity diagram of Fig. 3(a) 
and the tadpole diagram of Fig. 3(b). For the isospin propagator, one must 
sum over contributions from both pion and kaon intermediate states whereas 
the hypercharge propagator gets contributions from just kaons. 

It is the unitarity diagram which contains the 'real physics' such as the 
two-particle branch points. We shall explicitly consider the hypercharge case 
as it is somewhat less complicated. A straightforward calculation shows the 
invariant amplitude for the unitarity diagram of Fig. 3(a) to be 



One-loop Amplitudes 






unitarity 
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+ 6 iB2i{q^,m\) 



(19) 



The unitarity contribution cannot by itself be physical since it is noncovari- 
ant, divergent and scale-dependent. However, the tadpole amplitude 



tad 



?,iA{m\)g' 



(20) 



is seen to remove the problem of noncovariance in the unitarity amplitude 
via cancelation. The full hypercharge amplitude is then 



1— loop 



iA{m\) 
2mf^ 



6^52l((?^m^)+^: 



ct 



1— loop 



(21) 



where n|||i_ioop contains the counterterms generated by the next-order chi- 
ral lagrangian This latter contribution is shown in Fig. 3(c) and will 

be considered in Section 4. 

The above discussion has concerned determination of just the hyper- 
charge vector-current propagator at one-loop order. For completeness, we 
display without further discussion the one-loop amplitude for the isospin 
case (a = 6 = 3), 



n33('z' 



_ iA{ml) 
1— loop 3m^ 

,2 ™2 



4 iB2i{q\mi) - 2 ^S2l(g^ m^) + H 



1— loop 



(22) 



A simple check of the calculation is obtained by passing to the 5'f7(3) limit 
of = rriK- As expected, the isospin and hypercharge amplitudes agree. 

Two-loop Amplitudes 

In some sense, the content of this subsection must parallel that just 
given. One simply calculates the set of amplitudes associated with Feynman 
diagrams at two-loop order. However, in perturbative quantum field theory, 
things get worse as one proceeds to higher orders. Such is the case here. 
Instead of 2 types of non-counterterm amplitudes as in Fig. 3(a)-(b), there 
are now 9 non-counterterm amplitudes shown in Fig. 4(a)-(i). Of these, all 
but the 'eyeglass' diagram of Fig. 4(i) may be considered as corrections of 
one kind or another to the one- loop diagrams of Fig. 3. A detailed account 
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describing evaluation of these diagrams is clearly too lengthy a process to 
reproduce here. However, we can perhaps convey to the reader some useful 
insights with the following observations: 

1. Since the different vr, K and i] flavors can generally occur in each loop, 
there are many terms contributing to both the isospin and hypercharge 
two-loop amplitudes. 

2. A number of diagrams have noncovariant {g^u) contributions but they 
all turn out to cancel, as they must. 

3. Interestingly, some diagrams contain terms in which a divergent quan- 
tity is multiplied by a nonpolynomial function. Such terms are po- 
tentially disasterous to the calculational program since they cannot 
be removed by the usual renormalization procedure. However, they 
cancel totally. 

4. Two-loop amplitudes contain a number of contributions from the chiral 
lagrangian C^'^\ in addition to those from C^'^\ 

The final comment motivates introduction at this point of C^^^ , the next- 
to-lowest order chiral lagrangian. In principle, the entire 0{q^) lagrangian 
has three parts. The first has terms containing the meson field variable U, 

10 

= E ^iO^ 
1=1 

= Li [tV (^D^UD^'U^^f + L2TT [d^UD^U"^^ ■ Tr (d^'UD'^U^^ 

(23) 



£(4) 



+ L4 Tr (p^UD>'U^j Tr (xU^ + Ux^ 
+ L5 Tr (d^UD^V^ (xU^ + Ux^) ) + Lg [ Tr (xU^ + Ux^ 



Tr (x^U - Ux^)] + Ls Tr (xU^U^ + Ux^Ux^ 



+ iLg Tr [L^^D>'UD''U^ + R^^D^^UW^U 



The second has additional contributions occurring at the same order but 
which contain only external sources. 



Ce^t = Hi Tr {L^,L^"' + R^.R^") + H2 Tr (x^x) • 



(24) 
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In the above, Luy,Rny are the field-strength tensors 



L 
R 



(25) 



The third consists of two terms which vanish upon using the equations of 
motion, 



-Ceq-motn = GiTl {£■£) + G2Tt ■ i [x^U - x)) 



where 



£ ^ i (2D^ {U^D'^U) + x^U -U^X-\ Tr (x^U - U^x) 1 



(26) 



(27) 



In apphcations up to one-loop order, these terms can be neglected since 
they contribute only through tree diagrams where the equations of motion 
£" = may be applied. This is no longer true for a two-loop calculation. The 
terms in jCeq-motn may in general contribute a polynomial to the amplitude, 
which in turn could be absorbed into the coupling constants of the 0{q^) 
counterterm lagrangian.l&il In our application we have kept these terms 
at all stages of the calculation. Although present in intermediate steps, 
however, the effect of Gi, G2 cancels completely in the final result. 

For the hypercharge vector-current propagator the result of working out 
all the two-loop Feynman diagrams is 



n88('7' 



2-lp 
1 



iA{m\) 



QLio + (6 



m 



2 _2 



+ 3- 



q 

2 \\2 



m 



iB2i{q\mj^) + 18(^^21 ((?^m^)) 



K 



{ml + 2ml){U - 2Le) + 8m^(L5 - 2Ls 



-iA{m. 



+ n: 



2— loop 



(28) 



where we display only implicitly the two-loop counterterm contribution 
n8g|2-ioop of Fig. 4(j) and defer its calculation to Sect. 4. Observe how 
the above result reinforces the statement made earlier that all Feynman 
integrals in our analysis are reducible to just the two quantities A and S21. 
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To summarize, up to this point we have determined the isospin and 
hypercharge amphtudes through two-loop order, 112 



1— loop 



2— loop 



(a 



3,8) 



(29) 



These results cannot yet be accepted as physically meaningful due to the 
presence of divergent terms and to the explicit dependence on an arbitrary 
scale n which is an artifact of splitting loop and counterterm contributions. 



4 Enumeration of Counterterms 



The difficulties of divergence and scale-dependence in the 1-loop and 2-loop 
amplitudes are overcome by making appropriate use of counterterms. It 
is necessary to first generate a list at both 1-loop and 2-loop order of all 
counterterms which can possibly contribute. Carrying out this process for 
the 0{q^) counterterms will be the primary objective of this section. 

Let us begin by reminding ourselves of which counterterms appear at 
0{q^) order. a Out of a total of 12 possible operators in the chiral la- 
grangians C^^^ of Eq. (23) and £ext of Eq. (p^, only the Liq and Hi terms 



turn out to have nonzero matrix elements taken between single external 
vector-source lines. One easily finds the 0{q'^) counterterm amplitude to 
have the flavor-independent and covariant form. 



TTCt 

^^ab 



= -2{Lio + 2Hi) (a = 6 = 3,8) . (30) 

I— loop 



As noted earlier, the 0((?^) counterterm contribution appears in Fig. 3(c). 
Determination of the 0{q^) Counterterms 



Since the construction to follow will be somewhat formal in nature, it is 
important to keep in mind certain tenets of ChPT for guidance. Thus — we 
have found the two-loop calculations of 1133 and Ilgg to produce ultraviolet 
divergences. The low energy counting scheme of ChPT implies that these 
divergencies can be absorbed into the coupli ng c onstants of counterterms 
of order q^, thus renormalizing the amplitude. 113 However, more generally 
speaking all possible counterterms of order q^ have to be included in order 
to obtain the most general amplitude consistent with chiral symmetry. This 
is irrespective of whether divergencies occur in the loop amplitude or not. 
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In addition to those counterterms that contribute to 1133 and Ilgs, we 
shall also construct the set of 0{q^) counterterms which enter the two-loop 
analysis of the process 77 — > tt'^tt'^ and r] — > 7r''77. Our motivation stems in 
part from a recent analysis at two-loop order in ChPT of the first of these 
reactions. [[Lrl] In that work, however, no reference is made to a complete 
basis of 0{q^) counterterms. Our aim here will be to establish the precise 
relation between 0{q^) counterterms which contribute to both the 77 — > 
tt^'tt^ amplitudes and the vector current two-point function.^ In addition, 
we hope that the details given below will serve to simplify comparison with 
future works in this subject. 

Our first step is to derive a list of counterterms contributing a priori 
to the processes under consideration. Since we consider only couplings to 
neutral pseudoscalar mesons, the analysis simplifies considerably. There 
are no one-particle reducible diagrams involving 0{q^) counterterms. The 
reason is simply that to lowest order, q^, the vector source couples to an even 
number of Goldstone bosons. We may therefore switch off the external fields 
except for the vector sources v'^. Furthermore, since all flavour matrices 
involved are diagonal, the situation can be described by 

a^ = 0, p = 0, s = m, 
(f) = diagonal — > U = diagonal . (31) 

To proceed, we choose a set of building blocks transforming homogeneously 
under chiral transformations. The invariants can then be constructed easily 
by taking traces in flavor space. Our construction will closely follow the 



work in Refs. |lj, |15|]. The building blocks we use are 



Wfj.u = V^Uj. + S/uUf, , (32) 

as well as the totally antisymmetric tensor €^ypa- All the quantities U , x-, 
DjjU and the field strengths L^;^, R^j^^ have been defined in previous sections 
and u is the 'square-root' of the matrix-field U , 

U = u-u . (33) 



^Throughout this article we work in chiral SU(3). Comparison to the constants ap- 
pearing in Ref. p9[ is therefore only possible after the SU(2) limit has been taken. 
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Under a chiral SU{3) SU{3) transformation, the matrix u transforms as 



u{<P) ^ gnu{<P)h\^) = h{ct>)u{ct>)g[ , (34) 



where the nonhnear transformation h((t>) which first appears in Eq. (|34| ) is, in 
fact, defined by this same equation. ll^' 113 As a consequence, the left-chiral 
current transforms homogeneously, 

Uf,^h{^)u^{^)h^{(l)) . (35) 



Given some object O, a covariant derivative V^O which transforms in hke 

aed as 

Vf,0 = d^,0+[T^,0] , (36) 



manner to it^ can be defined as 



where 

Ti, = ^(^u^[d^-i{v^ + af,)]u + u[d^i-i{Vfj,-af,)]u^^ . (37) 

Here f^, are the matrix- valued external fields with spin-one introduced 
in Sect. 2. 

The building blocks of Eq. (32) as well as covariant derivatives thereof 
thus transform according to 

O ^ h{(t>)Oh^ {4>) . (38) 

The most general set of counterterms of order is then constructed by 
multiplying together the building blocks in Eq. (32) up to the appropriate 
order in the low energy expansion and then performing flavor traces like 

TriO.Oj...), Tr(OiOj. ..)Tr(Ofe. ..)... Tr(0/...) . (39) 

In Eq. (^), Lorentz indices are to be contracted and symmetrizations per- 
formed such that the discrete symmetries P and C are conserved. The 
transformation properties of the various building blocks under P and C can 
be found in Ref. We also note that all quantities in Eq. (32), except 

for X±) are traceless. 

The procedure just described for constructing the invariants of order 
leads to a set of terms which is in general redundant. In principle this set 
of terms can be reduced further by the use of the equations of motion, trace 
and epsilon identities 0:31 as well as by partial integrations. We shall not 
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elaborate further on such procedures, since the case we are interested in is 
much simpler and does not require that level of analysis. 

Let us return to the restricted situation as characterized in Eq. (|3lD. 
Under these circumstances the building blocks of Eq. (^) are transformed 
into 

X± = 25om(C/t ± U) , f+^, = 2{d^vy - d,v^) , 

= iu^d^u^ , f_^^ = . (40) 

w^lu = d^Uy + d^u^ , 

The only nontrivial objects at our disposal are therefore X+-, f+^^u find co- 
variant derivatives thereof as well as the totally antisymmetric e-tensor. The 
covariant derivative is also simplified, and is explicitly given as 

V^O = d^O-i[v^,0] . (41) 

Since only the diagonal isospin and hypercharge components a = 3, 8 are 
studied in this work, it follows that V^x+ ~ V^M = 0. In the same manner, 
we see that /+^;/ must be linear in the vector source v^, so two powers of 
fJr^lu are needed. As regards the operators X- and tf^jy, we note that X- 
involves only odd numbers of mesons and that if^jyis parity-odd and hence 
cannot lead to an operator of the form we are looking for. Finally, the 
remaining two powers of 'small' momenta as well as the two pseudoscalar 
meson fields can only be provided by two powers of n^j or one factor of x+- 
We are thus left with an expression of the form 



= ^2 (^iTr(VA/rvV+M.) + ^2Tr(V^/rvV+A.) 

+ TV (/+;../rx+) + Tr U+^..fT) Tr (x+) 
+ Kio Tr (/+^,/f Tx^up) + Tr (U^J'^^^'^Up) + ...] (.42) 



The ellipses in this equation denote all possible independent operators which 
can arise from the terms containing Kiq and Kn by changing the order of 
the building blocks f+^u and Up within the flavour trace and/or by taking 
two flavour traces. We shall refrain from writing down such contributions 
explicitly since they involve only diagonal flavour traces and thus lead to 
interactions of the type 

4e2 



^(6)' _ ptip 
-^0 



+ iK[29arg'"' + K[,g'^g'r) Tr (Q2) Tr {d'^cf>d-<p) (43) 
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In the above, we define the electromagnetic field strength 



F, 



fJ,U 



duA^ 



(44) 



where An is the photon field, and take 



eQAf, , Q = -diag(2,-l,-l) 



(45) 



It is now straightforward to determine the counterterm contributions to 
1133 and Ilgs which arise from the lagrangian £(6). We obtain 



^33 



and 



2— loop 



2— loop 



^0 



q\2Ki + K2) + 4.mlK3 + 4{2ml + ml)Ki (46) 



q^{2Ki + K2) + Am^K^ + 4(2m|^ + ml)K^ 



(47) 

where the lowest order expressions for squares of the meson masses (c/. 
Eq. (|lO|)) have been employed. Note that the contribution of counterterms 
to the difference, 1133 — Ilgs, is a constant depending only on K-^. Also, 
from the four terms Ki, ... K4 which may contribute a priori, only three 
independent combinations actually appear. The renormalization of coupling 
constants such as Ki, ... K4 as well as their properties following from the 
renormalization group equations will be discussed in Sect. 5. 

As for the two-meson two-photon processes, we first note a relation be- 
tween 77 — > Tr'^vr" and rj — > 7r''77. Upon evaluating the fiavour traces in 
Eqs. (p2|),(44) we easily see that as long as the pieces which violate the 
Zweig rule are neglected {i.e. those terms with two flavour traces), then the 
amplitudes for 77 vr'^vr^and r] — > 7r'^77 will depend on the same linear 
combination of the {i^jj.Ell Let us not pursue this line of reasoning any 
further, but instead concentrate in the following on the 77 
The 77 



Tr^vr" reaction. 



^0^0 ]2iatrix element can be written as 



(vr°(pi)7r0(p2)|7('?i,ei)7fe,e2))in = i(27r)V^)(P/-P.)e2e^^V (48) 



with 



ytMU 



A{s,t,u)T{"' + B{s,t,u)T^'' . (49) 

The tensors Tl'2 are defined in Ref. [jl9| , and we shall adopt the notation and 
conventions used there. It is straightforward to determine the contribution 
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of 0{q^) counterterms to the Lorentz invariant amplitudes A, B, and we 



This completes our discussion of the counterterm contribution to the 77 
7r°7r° reaction. Observe that the derivative lagrangian also contributes a 
constant piece to the amplitude A. The information contained in the two- 
loop expression for the vector propagator is thus not sufficient to fix the 
constant shift in A at this order of the low energy expansion. 

5 Renormalization Prescription 

Now that we have a complete set of 0{q'^) and 0{q^) counterterms, the 
next step is to formulate a procedure which allows the divergences to be 
subtracted and the scale-dependence removed from the calculation. 



There exist different ways of implementing these steps, e.g. minimal 
subtraction, etc. We have found it most convenient to subtract in the quan- 
tity A. This amounts to expanding the various counterterms in decreasing 
powers of A. In the following, let K denote any of the 0{q^) counterterms 
and L represent any of the 0{q^) counterterms. For the former, we write 



find 




(50) 



Subtraction Procedure 



n=2 



whereas the latter is given by 



—00 



L = lW(m)a' 



n=l 



= (/x) a + (/x) + (/x) a ' + . . .] . (52) 
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In our two- loop analysis, we shall actually need only the K^'^\ K^^\ K^^^ 
and the L^^\ L^'^\ L^~^^ counterterms. 

The definition of 0{q'^) counterterms used here has a simple association 
with that originally introduced by Gasser and Leutwyler in their papers on 
one-loop chiral perturbation theory. We restrict our attention to Ref. 
since, like this paper, it too employs an SU{3) basis. Then starting from 
Eq. (7.25) of Ref. [|, 

Lp-^^ = TiX + L^ (i = l,...,10) 

Ff-^) = AiX + H^ (i = 1,2) , (53) 

we see that our and {H^""^} are related to the Gasser-Leutwyler pa- 

rameters by 

= Ti and = (i = l,...,10) 

i^f^ = A, and Hf=Hl (i = l,2) . (54) 

Scale Dependence 

Having chosen a subtraction procedure, we must next determine the scale 
dependence of the counterterm coefficients. This is done by first noting that 
although the full counterterm is independent of scale, 

dK ^ dL 

jjL— — = and ^— = , (55) 
dfi dfi 

the {K^^^} and {L^")} must themselves generally be scale-dependent in or- 
der to compensate the scale-dependence of the prefactors, in Eq. (51) 
and in Eq. (|5^). Insertion of Eq. ( |5l] ) into the first equation of Eq. (55) 
yields an infinite sequence of differential relations for the 0{q^) countert- 
erms, 

12—— = 



dfi 167r2 



di^(o) 2 



dn IGvr^ 



(56) 



167r2 
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where, in terms of the Euler constant 7, we define 



C = -(ln(4^)-7 + l) 



(57) 



An analogous but distinct set of relations is obtained for the 0{q'^) coun- 
terterms by inserting Eq. (^) into the second equation of Eq. ( ^5| ) , 



djjL 



dfj, 

d L(-^) 
dfi 







167r2 



167r2 



c 

167r2 



(58) 



id) _ 



It is straightforward to sequentially solve relations of the type given 
above, and we find for the first few of the {K^""^}, 



(independent of fi) 



K(2) 

+ ^ In^ I ^ 



2 / /^o 



327r2 



whereas the sequence of {LW} be gins as 



L«(/.) 
lW(/.) 

L(-«(^) 



L«( 



(independent of /i) 



^(°)(.o) + |^ln 



+ ^ In^ I ^ 



2 / /^o 



647r2 



(59) 



(60) 
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Up to this point, the {K^"^} and {L^"'^} are arbitrary. In Sect. 6, we 
shall show how to constrain the counterterms which are relevant to our 
calculation by imposing conditions of finiteness and scale-independence. 

Effect of Mass Renormalization 

For the final topic in this Section, we consider renormalization of the 
meson mass parameters. This intermediate step turns out to considerably 
simplify the form of both the isospin and hypercharge vector-current prop- 
agators 1133 and Ilgs. For definiteness, we shall again explicitly study just 
llgg and thus restrict our attention to the renormalization of the kaon mass 
parameter. 

Recall in Ref. that for kaons, the relation between the bare mass 
parameter mx and the radiatively-corrected quantity Mk was shown to be 



K 



(61) 



where 



Ak = ^ 



7?2 
-^0 



'ml{Ml + 2Ml){Lf 



2L. 



(0)^ 



+8M^(4°^ -2L 




(62) 



For the two-loop amplitude, one can replace mj^ by everywhere since 
the error made occurs in still higher order. However, at one-loop order 
Eqs.(pl|),(|6^) must be utilized everywhere, i.e. both where mj^ appears 
explicitly and also where it occurs implicitly within a function. The effect 
of this substitution is 



Ussiq ,m^) 



1— loop 



n88('7',M, 



1— loop 



K) 



(63) 



This 



It can be shown that the additive term in Eq. ( pSj ) has the effect of cancelling 
the final term in the form for n88(Q^)|2-loop appearing in Eq. 
leaves us with the more compact expression. 



n88(r,Mi) 



2 -loop -^0 



+ 2M(M|; 



QLio + (6 



Mi 



+ q' 




+3 



iB2i{q',Mf,) + 18{iB2iiq^Mf,)) 



Ml 



(64) 
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The isospin amplitude behaves analogously, except now both pion and kaon 
mass renormalizations must be administered. 



6 Results 

We now have sufficient theoretical machinery to obtain our final representa- 
tions for the isospin and hypercharge vector-current propagators. The first 
step in this process is to expand the one-loop and two-loop expressions ob- 
tained earlier in powers of the singular quantity A. When combined with 
the counterterm contributions which were derived in the previous section, 
the full amplitude is rendered finite and free of the arbitrary scale /i by 
appropriate choices of the renormalization constants contained within the 
counterterms. The discussion becomes clearest if we first consider the one- 
loop case in some detail, and then touch on the high points of the more 
complicated two-loop sector. 



Renormalization at One-loop Order 



It is again convenient to refer to the simpler hypercharge one-loop am- 



plitude to illustrate the process. From Eq. of Sect. 3 and Eq. (|92| ) of 
App. A, we have 



n88(g' 



1— loop 



2Lio - 4/^1 

d-A 



2M2 



K) 



K 



+2L^^{Ml) + AHr{Ml] 



QiB2i{q^,Ml) (65) 



Note that we have adopted the convenient renormalization scale of for 
the counterterms L^iq and . Of course, any other choice is also allowed. 



but there will be compensating logarithmic factors as in Eq. (|60|). 

The next step is the crucial one. One removes simultaneously the diver- 
gent factor of A and dependence on the scale // with the choice 



24^+4^1 



(1) 



-1 



(66) 



This result is consistent with that Ref. B, where background field renor 



malization yielded the individual values L^q = — 1/4 and H\ 



r(l) 



-1/8. We 
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are then left with the finite and scale-independent one-loop result 



n88(g' 



= -2LSJ(Mi) - 4ijf (Mi) - 6^B2l(g^Mi) . (67) 

1— loop 

Renormalization at Two-loop Order 



Renormalization at two-loop order proceeds in like manner, but with 
many more terms. Since the underlying logic is the same, we shall not dwell 
on the many details of the analysis. However, it is of some interest to observe 
that the leading degree of divergence at two-loops is quadratic, 



n88((?^ 



2— loop 



32M|- 



(2) 



3 



(2) 



+ 



16Mj 



K. 



(2) 



+ 



(68) 



where we have exhibited only the most singular terms. The A divergences 
are removed by means of the scale-independent choices 



(2) 



and 



K. 



(2) 



K 



(2) 







(69) 



Analogous procedures for the remaining part of the amplitude turn out not 
to constrain the 0{q^) counterterms in terms of numerical values as in the 
above equation. Rather, relations are obtained which connect the 0{q^) 
counterterms with those from the O(g^) sector. 



4°)(^) + Li^o^(^) 



(0), 



(70) 



(0), 



This completes the list of constraints. 

Final Form of Vector-current Propagators 

We express our final result for 1188 in terms of the function iB2i and 
various counterterms, all fixed at scale M|-, 
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-2lSS(M 



AHf\Ml)-QiB2i{q\Ml 

-(0) 



+ 



^0^ 



18(ii?2i(g^M^))^ - 2AL'^>{Mf,)iB2i{q\Mf,)-P{Mf, 



Q{Ml) 



-^0 



2 
3^ 



(71) 



where we define 

P(^) ^ 



2K;°)(/.)+/^f (/i)+4-^)(/.) 



(72) 



The quantities P, Q, R are the three new independent counterterms intro- 
duced at two-loop order. Observe that Lg ^and L^q have been absorbed in 
these quantities. This is an example of a general property of the generating 
functional, as pointed out in Ref. |l9| . 

We can characterize the content of Ilgg as given above in terms of 'chiral 
counting'. Thus the first line on the right hand side is the one-loop result, 
whereas the remaining lines contain higher-order (two-loop) contributions 
which are proportional to the perturbative parameters of smallness q'^/F^, 
M^/Fq, M'^/Fq. For completeness, we also write down the analogous but 
more cumbersome result for 1133, 



^33(g^M2,Mi) 



1 



In 



Mi 



-2LSy(M2 



K) 



487r2 M2 



2\ o,- D 



+ 



mr{Mf,)-AiB2i{q\M^ 
iB2i{q\Ml)f+2{iB2i{q\Ml)f + 




+8^52l((?^M2)^52l((?^M|-) - P(Mi) + 



Lf(M2 



K> 



127r2 



In 



M 



K 



M2 



-{2iB2i{q\M^)+tB2i{q\Mf, 
AMI 



+ 



1 



247r2 Ml 



/?2 
-^0 



Q{Ml) + R{Ml) + ^^^^ 



K. 



+ Lf^{M' 



K) 



87r2 



In 



M 



K 



Ml 



(73) 



After some algebra, one finds that the expressions in Eq. (|7l|) and Eq. (|7i 
agree in the SU{2>) limit. 
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7 Applications 

The vector current two-point functions may be expressed in terms of the 
corresponding Kallen-Lehman spectral functions. Combined with the high 
energy behaviour which follows from perturbative QCD, the low energy rep- 
resentations of the two-point functions derived in the previous section is 
seen to contain two pieces of information: i) it yields the leading and next- 
to-leading order contributions to the low energy expansion of the spectral 
functions, and ii) it implies a set of sum rules which the spectral functions 
must satisfy. We shall discuss in turn the consequences of our results for i) 
and ii). 



An important application of our results involves determining the spectral 
functions associated with the isospin and hypercharge vector currents. As 
pointed out in the Introduction, the former is of special interest since it 
can be compared directly with existing data. Thus it is the isospin spectral 
function py33{q'^) that will be discussed here. 

From Eqs. @-(^, it follows that the relation between 1133(5^) and 
PV33(9^) is simply 



where sq = 4M^ is the threshold value. Thus the spectral function pyssiq'^) 
is proportional to the imaginary part of the propagator, 



Vector Spectral Functions 




(74) 



Im n33(g') = neiq^ - Ami) PV33{q^) ■ 



(75) 



A modest amount of effort in taking the imaginary part of Eq. (^) then 
yields the expression 
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-8Rei52l(g^ - 4Reifi2l(g^ M^) + 



(76) 



This result can be directly compared with data. Since such comparison 
is performed near the tttt threshold, only the part of the above formula 
proportional to 0{q^ — 4M^) contributes. Some care must be exercised in 
evaluating Re iB2i{q^ , M"^) and especially Re iB2i{q^ ,M]^). 

The one- loop result was already depicted in Fig. 1. A more complete 
comparison appears in Fig. 5, where we display the individual one-loop, 
two-loop, and total chiral predictions together with spectral function data. 
The two-loop result given in Eq. ( |76[ ) contains several additive contributions. 
Of these, the largest arises from the one containing the one-loop counterterm 
Lg'^^(M^). For our numerical work, we have adopted the value of L^^M"^) 
cited in Ref. |^], with a small correction for working at scale M^, 

= 0.0071 ± 0.0003 + 0.0002 = 0.0073 ± 0.0003 . (77) 

Another issue associated with a numerical analysis of /Ov33(q'^) is the value 
to be taken for the bare decay constant Fq (c/Eq. (|7^)). We have followed 
standard procedure by employing the physical value 

Fq ~ 0.093 GeV , (78) 

since any difference between Fq and F,r jWiH be relegated to still higher order. 
The effect of using the one- loop result o Fq < F^ would be to increase the 
two- loop contribution relative to that shown in Fig. 5. Finally, we caution 
that some care be taken in interpreting a graph such as the one in Fig. 5. For 
the purpose of illustration, we have depicted the energy-dependence up to 
the value s = 12 m^. However, before that value is reached one expects that 
still higher-loop effects should be taken into account. In addition, due to its 
stronger s dependence, the two-loop 'perturbation' ultimately becomes as 
large as the 'dominant' one-loop amplitude. From Fig. 5, it would appear 
that our calculation is valid up to s ~ 8m^ or an energy of about 400 MeV. 

Chiral sum rules 
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The real part of the vector current propagator depends on low energy 
constants appearing at 0{q^) in the chiral expansion. Some of these con- 
stants can be related to data by exploiting the high energy behaviour of 
QCD. We are thus able to derive a set of new chiral sum rules. 

The asymptotic behavior of the functions Hg niq '^) (a = 3, 8 not summed) 



is known from the operator product expansion, 



n..(52)~-^ 1 + ^^ In^ (a = 3,8), (79) 




where is the renormalization scale and as{q'^) is the QCD running coupling 
constant. It follows from the asymptotic form that the difference 1133 — Ilgg 
satisfies the unsubtracted dispersion relation, 

iU,,-Uss){q')= Tds (^V33- ms)(.) (3,) 
J so s-q^ -le 

The real part of this dispersion relation can be rewritten as a set of sum 
rules for the negative moments of the difference of spectral functions, 

d" (PV33 - pY88)is) I _ 
2v^(n33 - n88)(0) = n! / ds —-^ (n > 0) . (81) 



(dg 



so 



In this form, the relevance of calculating na(j(g2) at low is evident. 

Armed with the two-loop results obtained earlier in this paper, we can 
express the sum rules of Eq. ( pl| ) in terms of low energy constants of the 
effective chiral lagrangian. For instance, we have for the cases n = 0, 1, 

, (pV33-m8)(^) _ 1 

s ~ AStt^ 

AMI . Lr(Mi) +4g(Mi) ^ 16(Mi M^)^^^, ^^^^^ 



and 



87r2 M2 3F(2 



(ms - PV88)(s) 

ds 



so 



s 



2 



1 




18Fo2 Vl67r2y I M2 / Utt^F^ 



(83) 

The sum rule in Eq. ( |82D has already been obtained by Knecht, Moussallam 
and Stern. @1 However, since these authors calculate in generalized ChPT 
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up to 0{q^), only part of the symmetry breaking term proportional to Q ap- 
pears on the RHS of their sum rule, i.e. they necessarily lack the Lg and 
L^Q^^ dependence occurring the defining relation for Q in Eq. (^). By con- 



trast, our expression in Eq. (|82D is correct to 0{q ) in standard ChPT, and 



up to this order it is the combination of counterterms ^K^^ ~ ~ -^lo^'' 
which occurs. We note in passing that in addition to the term proportional 
to Lg'^ + -^10^ there would also be genuine two-loop contributions if the low 
energy constants were evaluated at a scale different from M^. Turning to 
the second sum rule, we see that it contains only one low energy constant, 
Lg*^\ Although l'^^ is known to rather high accuracy from the electromag- 
netic charge radius of the pion,lll it is nevertheless true that the relation in 
Eq. ( |83[ ) offers a nontrivial check on the determination of this quantity. 

Our two-loop representation of 1133 — Ilgs allows one to calculate the 
leading and next-to-leading order contributions to the RHS for all sum rules 
defined by Eq. (|8l| ) if n > 0. The leading term comes from the one-loop 
expression and is determined entirely by the vrvr and KK continuum terms. 
The next-to-leading order term gets contributions from both continuua and 
resonances, the latter being parametrized by low energy constants like Lg'\ 
It might well be that these contributions are numerically as important as the 
leading term since the relevant physical phenomenon enters the calculation 
first only at that level of approximation. 

Thus far, we have considered only the set of chiral sum rules associated 
with the difference, 1133 — Ilgs, of the polarization functions. We can also 
derive sum rules for the individual vacuum polarizations. According to the 



asymptotic behaviour given in Eq. (|79D , at least one subtraction is required 



to obtain convergent sum rules. Thus we can write 



00 



PYaais) 



ds ' ' = -— ^na„(0) (a = 3, 8) (84) 

for n > 1. As an example, we give the explicit expression for the isospin 
vacuum polarization a = 3 with n = 1, 

' PnM _ 1 / 1 1 \ 

s2 4807r2 Ml ^ 2Ml ' 



SF^ V16^V V 3 M2 y 87r2F(2 \^ 3 m2 



^P{Ml) . (85) 
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This sum rule clearly serves to fix the 0{q^) renormalization constant P. 
Note that the s"" sum rules which involve just the individual spectral func- 
tions are more slowly convergent than those involving a flavor difference of 
spectral functions because for the latter the leading term in the asymptotic 
expansion cancels. 

To analyze the sum rules of Eqs. (l8^)-(|85|), it will be necessary to eval- 
uate integrals containing negative moments of the spectral functions. We 
shall relegate such studies to a forthcoming publication, where we will also 
attempt to give realistic error bars on our determinations. Here, it is suffi- 
cient to touch upon the general strategy for calculating such integrals. As in 
Ref. 0, it is natural to divide the region of integration into three parts: A) 
the threshold region, B) an intermediate energy region where experimental 
information on the spectral functions is available, and C) the asymptotic 
region. The asymptotic behaviour of the spectral functions is calculated us- 
ing the OPE. As such, we anticipate the contribution from Region C to be 
especially small for the flavor-difference sum rules since for these the OPE 
will be proportional to light quark masses. Region B is the most difficult 
one to treat because it is here that intricacies of data analysis enter. Finally, 
in Region A we can use the improved low energy expansion of the spectral 
functions obtained in this paper. 

What will be the outcome of numerically analyzing the new chiral rules? 
The answer is that one thereby determines the combinations of counterterms, 
particularly those of 0{q^), appearing in them. By no means is this a trivial 
task, as the following example will demonstrate. For definiteness, we focus 
on the sum rule of Eq. (82) and rewrite it as 

{PY33 - PYS8)is) Lln.^^ 



so 



487r2 M2 

4M2 4°) (Mi) + 4? (Mi) ^ Ml , 16(Mi-M 



so that only terms proportional to the 0{q'^) and 0{q^) low energy con- 
stants appear on the RHS. We see from Eq. (^) that these counterterms 
are fixed by the difference between a spectral integral and a logarithmic 
mass ratio. However, the latter term can itself be expressed as a kind of 
spectral integral by employing the leading-ordei]^ chiral representations for 
the isospin and hypercharge spectral functions. To demonstrate this, let us 

^Similar arguments apply to next-to-leading order, but the corresponding expressions 
are much more complicated. 
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first define energy scales A3 and Ag below which the leading-order chiral rep- 
resentations Py23^ and Pygg^ are reasonable approximations to the physical 
spectral functions. We then carry out the integration to obtain 



A3 



ds 



AMI 



As 



ds 



487r2 



where the function /(sq/A) is given by 

3/2 



I 



so 



^ ds 
s 



1 



£0 

s 



2 


1 - 


So' 


3/2 

- 2 


1 - 


So 


1/2 


3 


A. 




A. 





AMI 
A, 



ln^ + 21n 
A 



1+1 



As 



(87) 



£0 
A 



1/2- 



Upon removing the cutoffs in an appropriate manner, one regains the loga- 



rithmic term appearing on the LHS in Eq. (86) 



VIZ. 



lim 

A3, 8-*°° 
A3-.Ag 



4M^ 
A, 



AMI 
As 



In 



M2 



(89) 



The point of this exercise is to stress that determination of the countert- 
erms will involve spectral integrals over A < s < 00 of the difference 
^^hys) _^(chir)^ rpj^^ reader should appreciate the irony that, although nega- 
tive moments enhance the threshold region in the new chiral sum rules, this 
has little influence on the counterterm determination. That is, as just shown 
the threshold region, in which p^^^^^ ~ Py^^^\ is suppressed via cancelation. 
Thus, to correctly infer values for the counterterms one must carefully ana- 
lyze regions in which p^^^^^ 7^ p^^^"^^ , most crucially the intermediate energy 
region where available data fix the spectral functions. 

Finally, let us point out that our discussion of Eq. (|8^) has yielded yet 
another insight. If one passes to the limit of, say, SU (2) chiral symmetry, 
then the explicit logarithmic dependence on in Eq. ( |86| ) leads to a di- 
vergence. In order to maintain consistency, then so too must the spectral 
integral in Eq. (^) diverge. The discussion just given shows precisely, via 
the chiral representations for the spectral functions, how this comes about. 
Thus the explicit logarithmic dependence in the chiral sum rules correctly 
reflects the leading implicit infrared singularities of the spectral integrals 
which appear therein. 
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8 Conclusions 



We have used ChPT to determine the two vector current propagators 1133 
and Ilgg through two-loop order. The main resuhs appear in Eq. (^TJ) and 
Eq. (|73|) . As expected, one finds the analytic structure corresponding to 
both TTTT and KK thresholds for 1133 and just the KK threshold for Ilgg. 
There is no Sir threshold for the propagators because we do not consider 
anomalous currents in this paper. The final results are seen to containa total 
of six renormalization constants and this warrants some discussion. The 
presence of the 0{q'^) counterterms L^q and already follows from the 
original one-loop analysis of Gasser and Leutwyler, whereas the remaining 
four arise from two-loop effects. Although the quantity Lg"^ is known from 
the one-loop analysis of the pion form factor, the other three (P, Q, R) 
are new. As we have seen, all but Lg'^^ contribute only to the real parts 
of the isospin and hypercharge propagators. The presence of these new 
renormalization constants is, of course, no surprise. For an effective theory 
like ChPT, the number of terms in the lagrangian increases sharply with 
loop order, and the renormalization constants are simply the set of quantities 
needed to render the theory meaningful at a given order. We shall comment 
shortly on how to partially constrain the new renormalization constants. 

There are a number of interesting applications of the work described here. 
Of these, the most immediate involves the comparison between data and the 
chiral prediction for the isospin spectral function pvsal^^)- Our calculation 
allows extraction of the isospin and hypercharge spectral functions at next- 
to-leading order in the chiral expansion. For /JvssC?^) the correction turns 
out to be substantial even at energies very close to the two-pion threshold, 
cf. Fig. 5. From this figure we can see that our two-loop representation of 
the isospin spectral function is valid up to an energy of about 400 MeV. 

Mention was made above of the counterterms which contribute to our 
final results. Three of these are new in the sense that they lie outside the 
phenomenology already established for the one-loop sector. They all occur 
in the real parts of the isospin and hypercharge propagators. By exploiting 
the asymptotic behaviour of the vector polarization functions 1133 and Ilgg as 
inferred from the operator product expansion, we have derived a set of new 
chiral sum rules. Among these is a sum rule for the 0{q'^) counterterm L^"^ . 
Since this quantity is already well-known, the sum rule offers a nontrivial 
check on the theory. The other sum rules fix 0{q^) counterterms, some 
of which also appear in the two-loop analysis of the process 77 tt^tt^. 
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We have begun the program of study outhned in Sect. 7 to constrain these 
quantities in terms of data, and results will be announced elsewhere. 

Finally, work is well underway on a two-loop determination of the SU{3) 
isospin and hypercharge axialvector current propagators. This study is a 
natural extension of both the present one and also of previous work done 
on the generalized Weinberg sum rules. M Although a chiral analysis of the 
axialvector propagators will superficially resemble that of the vector prop- 
agators, the two turn out to be very different in practice. In particular, a 
great deal of work on the renormalization of masses and decay constants 
of pseudoscalar mesons is required for the latter. Thus, we shall treat the 
axialvector sector separately in a forthcoming publication. 
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Appendix A: Feynman Integrals 

In this Appendix, we shall collect together definitions and properties of 
functions which play a central role in our analysis. In each case, one begins 
with an integral defined in d-dimensions and proceeds with its evaluation 
by isolating any divergence it might contain and determining the associated 
finite part. 

The Integral A 

The simplest such quantity is the scalar integral 

d'^k 1 



{2-kY k'^-w? 
The evaluation procedure is standard and we obtain 

A-d / ^ 



(4^)d/2 ^i-d 
d-A 
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(90) 

(91) 
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(93) 



The quantity /i, introduced in Eq. (piD, is the mass scale which enters the 
calculation via the use of dimensional regularization. The fi'^~^ prefactor in 
Eq. ( |92|) ensures that A{m'^) has the proper units in d-dimensions. 

The Integrals B, B^, B^^, 



Next we define 
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Since the scalar integral B is divergent, it is often convenient to work with 
a function B made finite by subtraction, 



B{q^,rr?) = B{q^, m^) - B{0, m^) 
where, with the aid of Eq. (p2|) , we have 



B{0,m'^ 



dm? w? 
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The subtracted integral can be written as 
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dx log 1 — xil — x) 







+ 



967r2 m? 9607r2 



+ ... 



(99) 



To study i?^, we invoke covariance to write B^{q^^m?) = -Big^. Then, by 
considering q^Bfj_ one finds Bi = B /2 so that 



B,iq\m') = '-^Biq^m') . 
The tensor integral Bn^ can be conveniently decomposed as 



B^t,{q^,-n?) = B2i{q^ ,w?)qf,q^ + B22{q^,m^)g^ 



(100) 
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For notational simplicity, we shall temporarily suppress function arguments 
in the following. Study of the contraction q^B^y leads to the relation 
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whereas the trace g^'^B^i, yields 

q^B2i +dB22=A + m^B 
From these relations, it follows that 
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A combination of these integrals which appears naturahy in loop dia- 
grams is 

T^u = ^B^y + q^qyB - 2q^B^ - 2qyB^ 

= AB^^-q^q^B (105) 
= {q^.qu - 9^uq''){^B2i - B) + 5^,(4^22 + ^q^B2i - q^B) , 

where we have introduced a convenient tensor decomposition for T^^ in the 
last line. The relations in Eq. ( |104D can then be used to obtain 



T^^u = {q^^qu - 9^.q^) (4i?2i + ^) + 2Ag^ 
where we define the finite quantity B21 as 
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The Integrals C, C^, C^,y 
Finally we consider the set of integrals 

C^u{q'^,m^) = 
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We first use covariance to write the vector integral as = Ciq^. Analysis 
of the contracted form q^C^ then determines Ci in terms of C, 



Ci{q\m') = ^ [q'C{q\m') + B{q\m^' 



(111) 



The scalar integral C can itself be determined by beginning with the identity 
q^dB/dq^ = —2q^C^, from which it follows 
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Analogous to our treatment of the tensor integral B^i, , we write the quantity 



C^i/ as 



As before, the contraction q^C^^ leads to one relation, 



q^C2i + C22 — ^ 
and the trace g'^'^Cuu yields another, 



q^C + B + B 



q^C2i + dC22 = B + m'^C . 



(113) 



(114) 



(115) 



Finally we define a quantity Unu which, like Tn^, occurs in loop diagrams, 
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With the input of the relations in Eqs. ( 1141) , ( pA^ ), the final line in Eq. ([IT 
results in the useful expression 
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Appendix B: The Isospin Vector Spectral Function 

In the following, we shall show how to obtain an expression for the vector- 
current spectral function good to two-loop order by using knowledge of the 
one- loop formulae for the pseudoscalar meson form factors. This will serve 
as a useful check of our two- loop analysis. For definiteness, we shall focus 
on the isospin spectral function. 

We begin with the statement of unitarity for the isospin vector-current 
propagator. At sufficiently low energies, only the physical two-pion inter- 
mediate state contributes. We can use Eq. (15) of Ref. and surrounding 
discussion to write 
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If we interpret both sides of this equation in the context of a chiral expansion, 
then the spectral function to (n+1)'^* order is expressible in terms of the pion 
form factor given to n*'* order. For example, in this way one immediately 
recovers the one- loop isospin spectral function of Ref. 
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The corresponding prediction for the spectral function at two-loop or- 
der is obtained in like manner, although requiring more work. The chiral 
expansion of the pion form factor appears to first order in Ref. |25], 

F^^{s) = 1 + 2H{s, Ml) + H{s, Mi) + . . . 

where, in terms of our notation, 
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It then follows to two-loop order that 
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Upon noting that 
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we obtain precisely the two-loop result of Eq. ( [76| ) derived earlier in this 
paper. 
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Figure Captions 



Fig 


. 1 


Comparison of 1-loop chiral prediction with data. 


Fig 


. 2 


Generic graph for computing vector-current propagators. 


Fig 


. 3 


One- loop graphs for Ay^^^. 


Fig 


. 4 


Two-loop graphs for ^b- 


Fig 


. 5 


Comparison of 2-loop chiral prediction with data. The one- 



two-loop and total contributions are denoted respectively by dotted, 
dashed and solid lines. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9501318v3 
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